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Abstract 

We estimate the power corrections (infrared renormalon contributions) to the coefficient 
functions for the transverse, longitudinal and asymmetric fragmentation functions in 
e + e~ annihilation, using a method based on the analysis of one-loop Feynman graphs 
containing a massive gluon. The leading corrections have the expected 1/Q 2 behaviour, 
but the gluonic coefficients of the longitudinal and transverse contributions separately 
have strong singularities at small x, which cancel in their sum. This leads to 1/Q 
corrections to the longitudinal and transverse parts of the annihilation cross section, 
which cancel in the total cross section. 
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1 Introduction 



Experimental studies of the single-hadron inclusive spectrum in the e + e~ annihilation process 
e + e~ — > 7*/Z° — > hX have been performed with high precision over a wide range of energies 
for various types of produced hadrons h (see Refs. and references therein). These studies have 
mostly concerned the total fragmentation function, 

1 <!<T >e+e~ ^hX) = F tot (x,Q 2 ) , (1.1) 



(Ttot dx 



where x = 2Q-ph/Q 2 is the energy fraction of the observed hadron and Q M is the virtual boson four- 
momentum, equal to the overall centre-of-mass momentum in this process. The most precise data are 
those for the total fragmentation function into charged hadrons. The scaling violations (logarithmic 
Q 2 -dependence) in this quantity predicted by perturbative QCD |5|-|| have been observed and used 
to measure the strong coupling constant [@||. 

The ALEPH || and OPAL Q collaborations at LEP have also presented results on the joint 
distribution in the energy fraction x and the angle 9 between the observed hadron h and the 
incoming electron beam. We can write |3| 

' |(1 + cos 2 9) F T (x, Q 2 ) + \ sm 2 6F L (x, Q 2 ) + \ cos6F A (x, Q 2 ) , (1.2) 



(Ttot dx d cos #8 4 4 

where Ft, Fl and Fa are respectively the transverse, longitudinal and asymmetric fragmentation 
functions. As their names imply, Ft and Fl represent the contributions from virtual bosons polar- 
ized transversely or longitudinally with respect to the direction of motion of the observed hadron. 
Fa is a parity-violating contribution which comes from the interference between the Z° and pho- 
ton contributions. Integrating over all angles, we obtain the total fragmentation function (|1.1|) , 
Ftot = Fp + Fl- 

According to the factorization theorems of QCD ||, each of the functions Fp (P = tot,T,L or 
A) can be represented as a convolution of universal parton fragmentation functions Di (i = q,q or 
g) with perturbatively calculable coefficient functions C P 



F P (x, g 2 ) = E / — (*) Dfa/z, Q 2 ) . (1.3) 

■ Jx Z 

The parton fragmentation functions themselves cannot be computed perturbatively, although their 
logarithmic Q 2 -dependence, which is the main source of scaling violation, is predicted by the QCD 
evolution equations || . The evolution kernels (splitting functions) differ in non-leading orders ||[7| 
from those for deep inelastic structure functions. Thus scaling violation in fragmentation provides 
an important independent test of QCD. The longitudinal and transverse fragmentation functions 
can be used to measure the gluon fragmentation function while the asymmetric part may 

be useful for the measurement of electroweak couplings |Hj]. Once measured and parametrized 



PJ8] JTI|JT2] ] , the parton fragmentation functions can also be used to test QCD in jet fragmentation 
in other processes such as lepton-hadron Jl3| and hadron-hadron collisions. 



In addition to the logarithmic Q 2 -dependence predicted by perturbative QCD, it is expected 
that fragmentation functions will exhibit process-dependent power corrections, i.e. contributions 
proportional to inverse powers of Q, analogous to the higher- twist contributions found in deep 
inelastic scattering ||15|| . An understanding of power corrections is crucial for precision tests of 



i 



scaling violation. In the deep inelastic case such contributions are related to the hadronic matrix 
elements of local operators, and it is well established that the leading corrections should be of order 
1/Q 2 - The machinery of the local operator product expansion is not applicable to fragmentation, 
but a study of the relevant non-local operator matrix elements again suggests a 1/Q 2 behaviour 
|16| |. The current data on the total e + e~ fragmentation function F tot are consistent with either a 
1/Q or 1/Q 2 form for the leading power correction j3|. Data on the separate functions Fp for P=T, 
L and A are at present limited to a single energy Q = Mz HU , and so predictions concerning their 
scaling violation and power corrections remain to be tested. 

In the present paper we estimate the power corrections to fragmentation functions using a 
recently developed 'dispersive' method based on the infrared properties of Feynman graphs and 



some assumptions about the strong coupling at low scales [[lTj] . The techniques involved are similar 
to those involved in the study of infrared renormalons [|T^-|2T]|, which correspond to unsummable 
divergences of the perturbative expansion. Here we use them as a more general probe of the influence 
of soft regions of integration on hard process observables. The results obtained by applying these 
techniques to deep inelastic scattering |T7|j2^j23[| are consistent with those of the operator product 
expansion and look promising phenomenologically. 

A calculation of the leading power correction to the total fragmentation function F to t using the 
dispersive approach was presented in Ref. [17]. A 1/Q 2 correction was found and the corresponding 



quark coefficient function was computed. In the present paper we extend these results to subleading 
(1/Q 4 ) power corrections, and to the transverse, longitudinal and asymmetric quark fragmentation 
functions separately. We also compute (in a certain approximation) the corresponding gluonic 
coefficient functions. 

We shall be particularly concerned to clarify an apparent paradox which arises when one consid- 
ers sum rules for fragmentation. Summed over all particle types, the total fragmentation function 
satisfies the energy sum rule, which we may write as 

1 fl dxxF tot (x,Q 2 ) = 1 . (1.4) 



2 Jo 

Similarly the integrals 



1 C dxxF^{xM 2 ) = — (1-5) 

2 Jo otot 

give the transverse and longitudinal fractions of the total cross section. The perturbative prediction 



is H| 



ol ctt o(s /601 6 37 \ fa 



Otot CT t ot TT 



(1.6) 



~ — + (13.583 - 1.028n f ) ( — 

7T V 7T 

where £(3) = 1.202 and rif is the number of active flavours. Note that the whole of the 0(a s ) correc- 
tion to cr tot comes from the longitudinal part, while the 0(a 2 ) correction receives both longitudinal 
and transverse contributions. 

The OPAL data point || for ai,/a tot is shown, together with the perturbative predictions, in 
Fig. |I[ The data lie somewhat above the next-to-leading-order prediction (dashed), which suggests 
that higher-order and/or non-perturbative corrections are significant. An estimate of the latter is 
provided by the difference between the JETSET Monte Carlo |25j hadron-level prediction (dot- 
dashed) and the perturbative result. This difference shows a clear 1/Q behaviour, with a coefficient 
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Figure 1: Longitudinal fraction of the e + e hadronic cross section. 



of about 1 GeV. Assuming that the same behaviour will be manifest in the data, we have an apparent 
paradox: the function xF^ has only a 1/Q 2 correction, but its integral has a 1/Q correction. 

Power corrections proportional to 1/Q, usually referred to as hadronization corrections, are 
typical of hadronic event shapes in e + e~ annihilation, where they are clearly seen in the data and 
are predicted by a variety of approaches p6|-|2"8|j. String-like models of hadronization, for example, 
in which the energy of a jet is redistributed with proper density A ~ 500 MeV per unit rapidity and 
limited transverse momentum relative to the jet axis, imply a correction to the longitudinal cross 
section of [||[] 

Sa L ttA 0.8 GeV 

= — ~ . (1.7) 

Otot 2Q Q 

Adding this to the next-to-leading-order prediction gives the solid curve in Fig. Ill which agrees 
well with the JETSET prediction. The correction arises from mixing between the transverse and 
longitudinal angular dependences in Eq. ( |1.2| ) due to hadronization. The transverse cross section 
receives an equal and opposite correction, and so there is no 1/Q term in the total cross section. 



We shall see that dispersive approach of Ref. |17fl leads to the following resolution of the paradox. 
The 1/Q terms arise from soft gluon fragmentation. The gluonic coefficient functions of the 1 / Q 2p 
power corrections are highly singular at small x. Upon integration they are all 'promoted' to a 1/Q 
behaviour and have to be resummed. The result is a 1/Q correction to ol /ototj with a coefficient 
similar to that in the hadronization models. 

In the remainder of the paper, we first give a brief summary of the relevant assumptions and 
results from Ref. |17|. The dispersive method is based on the evaluation of one-loop Feynman 



graphs containing a gluon of finite mass fx. This yields the 'characteristic functions' for the relevant 



^Note that there is a misprint in the corresponding equation (3.35) of Ref. ||] 



quantities, which are given in Sect. 3. We extract the power corrections from the behaviour of the 
characteristic functions as p —>■ 0. The rules for taking this limit are also explained in Sect. 3. 
In Sect. 4 we give the expressions thus obtained for the power corrections to the fragmentation 
functions, and to the transverse and longitudinal cross sections. Finally in Sect. 5 we discuss the 
results and give some numerical predictions. 



2 Dispersive method 

We do not repeat the discussion of Ref. [|17|] but simply summarize the results required here. The 
basic assumption is that the dominant non-perturbative contributions to the coefficient functions 
Cp in Eq. (|1.3|) are of the form 

5Ci>(x, Q 2 ) = r % 5a eS (p 2 )^(x, Q 2 ] t?) • (2.1) 

J fi Z 

Here 5a e s(p 2 ) is a non-perturbative modification to the effective strong coupling, restricted to 
the region of low p 2 . T P (x,Q 2 ; p 2 ) is the relevant characteristic function, that is, the parton 
fragmentation function computed at order a s using a non-zero value p for the gluon mass in the 
Feynman denominators of the contributing graphs. T P is minus the logarithmic derivative of T v 
with respect to p 2 . Since JF P depends only on dimensionless ratios, we may write 



p. 2 

J-p(x,g 2 ;/i 2 )=^ P (x,e), ^(x,e) = -e-r p (x,e) , ^ = ^. (2.2) 



The non-perturbative contributions SCp thus depend on the small-e behaviour of JF P , which is of 
the generic form 

J*(x,e) = -P£(x)\ne + Cl P (x) - Q P (x)elne - |Q iP (x)e 2 lne , (2.3) 

where the dots indicate terms that are either (9(e 3 lne) or analytic and vanishing at e = 0. Here 
Pp(x) (contributing for P=T,A only) is the q — > % splitting function, Cg P (x) is the relevant per- 
turbative coefficient function (in the gluon mass regularization scheme), and G\ v {x) etc. will be 
related to non-perturbative corrections. 

A crucial point is that, for consistency with the operator product expansion, the integer /i 2 - 
moments of the coupling modification should vanish: 

^fPSaaQi 2 ) = , (2.4) 
/i 2 

at least for the first few moments p = 1, . . . ,p max ~ 9. As a consequence, only those terms in the 
small-e behaviour of J-p(x, e) that are non-analytic at e = lead to non-perturbative contributions 
19| . Thus from the small-e behaviour (|2.3| ) we find 



5CUx,Q 2 )=Cl P (x)^ + Q, P (x)^f + ••• , (2.5) 



where, following Ref. ||17|| , we have defined the log-moment integrals 

C F f°° dp, 2 



A' 2p = ^ jH ^ ^ ln(/i 2 //i 2 ) 6a eS ^ 2 ) . (2.6) 



i 



Notice that since integer /z 2 -moments of Sa e s vanish, these quantities are independent of the scale 
Hq. For convenience, we extract a universal factor of Cp/2-K from the characteristic function. 

Instead of interpreting the coefficients A' 2p in terms of a universal low-energy effective coupling, 
one may treat them more generally as pro cess- dependent parameters to be determined experimen- 
tally. Eq. fl2.5p still has predictive power because the coefficient functions C\-p(x) etc. specify the 
x-dependence of the power corrections. This dependence is supposed to reflect the relative sensi- 
tivity of different regions of x to soft dynamics. 

The technique we use to evaluate the coefficient functions of power corrections, viz. extraction 
of the non-analytic terms in the massive-gluon expressions for observables as fi 2 —>■ 0, is the same 
as that applied in studies of infrared renormalons In the language of renormalons, the terms 
computed are ambiguities in the perturbative prediction for the observable in question, which have 
to cancel against corresponding ambiguities in power-suppressed non-perturbative contributions. 
Here we argue that the non-perturbative contributions themselves should display the same power 
behaviour and x-dependence, since the small-// 2 limit probes the sensitivity of an observable to the 
soft non-perturbative region as a function of x and Q 2 . 

We expect quark and gluon fragmentation to contribute to power corrections on an equal basis, 
since the dominant contributions are assumed to be determined at first order in a e g. The application 
of the dispersive method to compute the quark contribution is straightforward since in that case we 
sum inclusively over all gluon fragmentation products. The inclusive sum generates a contribution 
which is equivalent to that of a massive gluon, as discussed in ref. ||17|| . We shall also use the 
dispersive approach to calculate gluonic contributions, but its application there is more questionable. 
By definition we observe the fragmentation products of the gluon in that case, which spoils the 
equivalence to a massive gluon. The situation becomes similar to that for an event shape variable: 
the correction has the same leading power behaviour as that due to a massive gluon, but the 
coefficient may be modified. This can be investigated in the large-nj limit, in which the gluon 
fragments only into quark-antiquark pairs |28| , |29| . Although the question requires further study, we 
assume here that the massive-gluon technique does provide a reasonable estimate of contributions 
from gluon as well as quark fragmentation. 



3 Characteristic functions 

The object of central importance in the dispersive method is the characteristic function J-p(x, e) 
for the emission of a gluon with mass-squared fj 2 = eQ 2 at the hard scale Q 2 . This is computed 
from the relevant one-loop graphs with a modified gluon propagator. The characteristic function 
for the total fragmentation function is obtained by contracting the resulting hadronic tensor with 
the tensor representing a sum over virtual-boson polarization states, 

O^O v 

E^p = -^ + ^- (3-1) 
p v 

The corresponding tensor for the longitudinal part is where e£ is the polarization vector along 
the direction of p h , the three-momentum of the observed hadron in the virtual-boson rest frame: 



The characteristic function for the transverse fragmentation function is then obtained as the differ- 
ence between the total and the longitudinal part. For the asymmetric fragmentation function we 
take the vector-axial interference term and contract with the tensor 

2 W \p h \Q 



\^HVcrp | I /"i ' (3-3) 



The resulting expressions for J-"p(x, e) are given below. The power corrections are deduced from the 
non-analytic terms in the small-e behaviour of the logarithmic derivative J- p (x,e). When taking 
the small-e limit, one should be careful with the phase-space boundaries, and in particular with 
functions that are singular on or near these boundaries. We give the rules for obtaining the correct 
limiting behaviour. 



3.1 Quark fragmentation 



The phase space for fragmentation into a quark of negligible mass with emission of a gluon of 
mass-squared eQ 2 is < x < 1 — e. The characteristic function is therefore of the form 

J%(x, e) = F£\x, e) 0(1 - x - e) + ^ v \e) 5(1 - x) (3.4) 

where J-^ ' represents the virtual contribution. For transverse quark fragmentation the real gluon 
emission part is 



(x + e)(l — x) 



2 1 + e 2 e e z , 

! --l-x-2e + A- + Q— In 

1 — X XX 2 

3-5e 2 e e 2 _e e(l - e) I 

2(1 -x) + (1-x) 2 + 2(1 -x) 3 + x + e x 2 ( x > + 6e - 



(3.5) 



The virtual contribution is 

T {v) {e) = 2(l + e) 2 
where 



1 7T 2 

Li 2 (-e) + In e ln(l + e) - - In 2 e + — - 

2 o 



Li 2 (-e) 



-In (1 + t . 

o t 



(3 + 2e)lne-2e, (3.6) 



(3.7) 



The characteristic function for longitudinal quark fragmentation has only a contribution from 
real gluon emission, which is 



Fl(x,e) = -2- 2 + 3- In 



(x + e)(l — x) 



e(l + 2e) e 2 e(l - e) 
2 \ ; + + 6^ '- + 1 - 2e . 



(3.8) 



1 — x 



1-x) 
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The characteristic function for asymmetric quark fragmentation is of the form 
same virtual contribution, but the real gluon emission part becomes 

"2(1 + e) 2 



1 — x 

3(1 ~6 2 ) 
2(1 -x) 



1 - x - 2e + 2 



c(2 + e) 



x 



In 



e e 
+ 



1-x) 2 2(1 -x) 3 x + e 2 



(x + e)(l — x) 



(1 -x) + e 



with the 



(3.9) 



In taking the small-e limits of Eqs. ( |3.5| ), ( |3.8| ) and Q3.9D , we must remember that the phase-space 
boundary is at x = 1 — e. Now for any function F(x) that is analytic in a neighbourhood of x — 1 
and any test function f(x), we have 

r e F(x)f(x)dx= f 1 F(x)f(x)dx- eF(l)f '(1) + \e 2 [F'{l)} '(1) + F(l) /' '(1)] + ■ • • • (3.10) 
Recalling that 

5 (n) (l - x) f{x) dx = f (n \l) , (3.11) 



we can make the e-dependence explicit, up to terms of order e 2 , by replacing any expression F(x) 
analytic at x = 1 as follows: 

F(x) -> F(x) - e[F(l) - §eF'(l)] 5(1 - x) + ±e 2 F(l) - x) . (3.12) 

For expressions that are singular at x = 1, we define '+', '++' and '+++' prescriptions such 
that, for any test function f(x), 

1 F(x)+ f{x) dx = f 1 F(x) [f(x) - f(l)] dx 
Jo 



F(x) ++ f(x) dx = / F(x) [f{x) - f(l) + (1 - x)f'(l)] dx (3.13) 

J 

F(x) +++ f(x) dx = £ Fix) [fix) - /(l) + (1 - x)fil) - |(1 - x) 2 f (1)] rfx . 
Using Eq. (|3.11 ), we can now replace the singular terms, up to terms of order e 2 , as follows:| 

\ne6(l-x) + e6'(l-x) - \t 2 8"{\-x) 
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il-x] 
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J - \ In 2 e 5(1 - x) + e(lne - 1) 5'(1 - x) - ie 2 (lne - |) <5"(1 - x) 
+ (1 - e)<J(l - x) + e In e <5'(1 - x) - \e 2 8"{1 - x) 
- + |(1 - e 2 ) 5(1 - x) - e(l - e) <5'(1 - x) - |e 2 lne5"(l - x) . 

(3.14) 

Applying these rules to Eqs. (|3.5|) etc., we obtain expressions of the form ( [2.3|) . The coefficient of 
— In e in A ^ s ^ ne q uar k splitting function P qq ix) = (1 + x 2 )/(l — x), which is singular for x — ► 1. 
The singularity is regularized by including the virtual contribution. As explained in Ref. flr7| , this 
term produces the logarithmic scaling violation in the quark fragmentation function. The second 
term Cq P (x) is the perturbative coefficient function (in the gluon mass regularization scheme). The 
remaining terms generate power corrections of the form (|2.5|), which we list in Sect. 4. 



3.2 Gluon fragmentation 

The characteristic functions for gluon fragmentation depend on the quantity p = \J x 2 — 4e and are 
defined in the phase-space region 2-y/e < x < 1 + e. Naturally, they have only real gluon emission 

'•Note that these rules are different from those given in Ref. j2^] for the deep inelastic case, since the phase space 
is different. 



contributions. For transverse gluon fragmentation we find 
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and for longitudinal gluon fragmentation 
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(3.15) 



(3.16) 



There is no gluonic contribution to the asymmetric fragmentation function. 

Since the upper phase-space boundary is now at x — 1 + e, the rule ( |3.12| ) for obtaining the 
small-e behaviour of expressions that are analytic at x = 1 becomes 



F(x) -> + e[F(l) + f eF'(l)] 5(1 - x) + \e 2 F(l) 5\1 



x) 



(3.17) 



There are no gluonic contributions that are singular at i = 1. Instead, the lower phase-space 
boundary, x = 2^/e, is e-dependent, and there are terms that are singular at or near this boundary. 
However, for any finite x the region of integration in Eq. (|1.3|) does not extend to the lower phase- 
space boundary, and so it and the nearby singularities are irrelevant. Thus for any finite x we can 
safely expand p = x — 2e/x — ■ ■ ■ and use Eq. (|3.17|) to obtain the small-e limits of Eqs. ( [3 .15 ) and 

When taking moments of the fragmentation functions, on the other hand, we integrate all 
the way down to x — 0. For sufficiently high moments the contribution of the small-x region is 
suppressed and the above procedure will still be reliable. For lower moments, the x-integration 
must be performed first, and then the small-e limit can be taken. We shall see that the phase-space 
boundary and singularities at small x can play a crucial role in this case. 



4 Power corrections 
4.1 Quark fragmentation 



The coefficients in Eq. ( [2.5|) for the first two power corrections to the transverse quark coefficient 
function are found from Eqs. (|3.4j )- (|3.6| ) to be 



(4.1) 



CJtO*) - (1 _ x) 



' -2 + - + 26(l-x)-5'(l-x) 

x 



Cj T (x) = + - 2 + 5 5(1 - x) + \ 5"(1 - x) 

1 — x i x z 2 



The corresponding expressions in moment space, defined by 



C(N) = / x"- L C(x)dx 
Jo 



(4.2) 



are 

q T (iV) = -Ar + 3-| + ^_-4S 1 



with 



CI, T (N) = ~iV 2 - \n + 6 + - 45! , 



N-l i 



For the longitudinal quark contribution, the corresponding results are 

Cl L (x) = ---2 5(l-x) 

x 

Cl L (x) = ~-86{l-x)- 26' (1 - x) , 



x' 2 



CI AN) 



N - 1 

Cl L (N) = -2N-6 12 



N-2 

For the asymmetric coefficient function, which receives only a quark contribution, 
(%A X ) = /-, 4 \ - 2 + - + 25(1 - x) - 5\l - x) 

1 — X L. X 



C! )A (x) = + ^ + 35(1"^) + ^ - x) 



C{ A (N) = ~N + 3 - | + - 4& 



(4.3) 



Si = V - = V(^)+7E = lnA^ + 0(l/iV). (4.4) 



(4.5) 



(4.6) 



(4.7) 



(4.8) 

C1 A (N) = lN 2 -^N + 4+ - 45! . 

Thus the 1/Q 2 corrections to the transverse quark and asymmetric coefficient functions are the 
same, but the 1/Q A corrections (and higher power corrections) are slightly different. 

Note that the expressions given above for the moment coefficients C\ vV with P=T,L are only 
correct for N > p. As discussed earlier, for lower moments the low-x singularities of the character- 
istic functions have to be taken into account, and the singularity structure in e becomes different. 
An important case is that of the N = 2 moments, which define the contributions to the transverse 
and longitudinal cross sections. This will be examined more fully in Sect. 4.3. 



4.2 Gluon fragmentation 



The coefficients in Eq. ( [2.5|) for the first two power corrections to the transverse gluon coefficient 
function are found from Eq. ( p. 15 ) to be 
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For longitudinal gluon fragmentation, the corresponding results are 
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(4.10) 



(4.11) 



(4.12) 



The above expressions for the gluonic moment coefficients C^ p p are only valid for N > 2p + 1. 
For gluonic moments with iV < 2p + 1 the phase-space boundary and low-x singularities of the 
characteristic functions become important and change the singularity structure in e, as will be 
illustrated in detail for N = 2 below. 



4.3 Transverse and longitudinal cross sections 

Defining the characteristic functions for the various contributions to the cross section as 7Zp(e) for 
P=T,L and i — q,q, g, such that 
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(4.13) 



where <To is the Born cross section, we would expect from Eq. 
corrections would take the form 



that the corresponding power 
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(4.14) 



where the coefficients C2 P (2) and C*4 P (2) are the N = 2 moment coefficients given above. This 
works more or less as expected for the quark contributions, yielding 1/Q 2p power corrections, modulo 
logarithms. The expressions given for Cf P (iV) in Eqs. ( f4.3| ) and ( (4.6| ) have a pole at iV = 2; this 
singularity corresponds to the logarithmic divergence of the x-moments of the expressions for Cf P (x) 
in Eqs. (|4.1|) and (4J3). Performing the x- integrations first, and then taking the small-e limit, we 
find that the quark characteristic functions are in fact 
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2 97 1 

- hie H h 3elne + 3e 2 lne + -e 3 lne + • 

3 36 3 



(4.15) 



i n 



In Eqs. (|4.15|) - (|4.18|) , the dots represent terms that vanish as e — > 0, and are either analytic or 
(9(e 4 lne) at e = 0. Note that 1Z\ has a lne divergence, since the separation of and a\ is not 
collinear safe. 

The (9(elne) term in gives a leading power correction of order 1/Q 2 in of. The effect of the 
logarithmic divergence of the x-moment of C\-p{x) is to 'promote' the (9(elne) terms to 0(eln 2 e), 
resulting in a InQ 2 enhancement of the 1/Q 4 corrections to and of. 

In the case of the gluon contributions, the relevant coefficients, given by Eqs. ft4.9| ) and (|4.11| ), 
are so singular as x — > that the 1/Q 2p power corrections are promoted to 1/Q. In the case of of, 
for example, we have 



,,2 „l +/ ,2/Q2 o „l J 



Similarly, the most singular parts of C| L (x) and all the higher coefficients give contributions of 
order 1/Q. Hence these terms have to be resummed, and the true behaviour of the characteristic 
functions for a^ T involves y/e singularities. Again performing the x-integrations first, and then 
taking the small-e limit, we find 

4 44 / 8\ 2 

TZ 9 T (e) = --lne- — + vr 2 (l + e) 2 ^ - elne(ln e + 4) - e 2 In e (lne + - +— e 3 lne + --- 
3 9 \ 6/ 15 

/ 1f(\ 92 

n 9 L {e) = 1 - tt 2 (1 + e) 2 v^+ elne(lne - 2) + e 2 lne ( In e - ^ ) - ^e 3 In e + 

n 9 tot (e) = -- In e - — - 6e In e - 6e 2 In e - -e 3 In r + • 
3 9 3 



3 J 15 



(4.17) 

Notice that the \fi singularities cancel in the total gluonic contribution. Adding the quark and 
gluon contributions together gives 

TZ T (e) = 7T 2 (l + e) 2 v / e-elne(lne + 4) + 2e 2 lneflne-— ) + -e 3 lne + --- , 

V 6/5 

3 / 17\ 22 

K L (e) = - - tt 2 (1 + e) 2 ^ + e In e(ln e + 4) - 2e 2 In e In e - — ) - — e 3 In e + • • • , (4.18) 

2 V d / 15 

^tot(e) = ^-^ 3 lne + --- , 



corresponding to 



( 3 ^ a s -n 2 A{\ . 
^--»(^ F - + T -) (4.19) 



where j4' 6 is defined by Eq. (|2.6j) with p = 3 and 



C F f°° dp 2 2 C 



Ai = — / — // 5a eff (/i ) = — / dp, 5a efi (p ) . (4.20) 
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X 

Figure 2: Coefficient of 1/Q 2 correction to the total fragmentation function (solid), with quark 
(dashed) and gluon (dot-dashed) contributions shown separately. 



5 Discussion 

To illustrate the above results, we have computed the predicted 1/Q 2 contributions to the total, 
transverse and longitudinal fragmentation functions using the parametrizations provided by the 
ALEPH collaboration |3[] for the parton fragmentation functions Di(x,Q 2 ). We express the results 
in terms of the coefficients -D 2 ,p where 



F P (x,Q 2 )^Fr\x,Q 2 )[l 



D 2 ,p(x,Q 



2\ 



Q 



2 



A' r 1 rl? 

D *A*, Q 2 ) = FPert , 2 n2 , E / -C\,A Z ) AO*/*, Q 2 ) 

Pt> (X, U ) „• Jx z 



(5.1) 



(x,Q 



the coefficient functions C\ v being as given above. In each case, Fp eit (x,Q 2 ) represents the rel- 
evant lowest-order perturbative prediction. We assumed the value A' 2 = —0.2 GeV 2 for the non- 
perturbative parameter defined by Eq. (|2.6|) with p = 1. This value is suggested by deep inelastic 
data P^| . Calculations were performed at Q = 22 GeV, but the Q-dependence of the coefficients 
D2,p(x,Q 2 ) is small. 

Figs 0, |3| and f| show the results for F tot , F T and F L , respectively. Note that the longitudinal 
coefficient is much larger than the others simply because we divide by Fl eIt , which is of 
order a s relative to F^ and F^ er . We see that, as discussed above, the contributions from gluon 
fragmentation give large opposing corrections to Ft and F^ at small x, which tend to cancel in F tot . 

Turning to the results on the transverse and longitudinal cross sections, Eq. fl4.19|) together 



with the empirical value (|1 . 7|) for the 1/Q correction to <7l would suggest a value A\ ~ 0.2 GeV for 



1.5 




Figure 3: Coefficient of l/Q 2 correction to the transverse fragmentation function (solid), with quark 
(dashed) and gluon (dot-dashed) contributions shown separately. 




1 9 



the non-perturbative parameter ( f4.20| ). 
deduced event shape data in Refs. 
was given as 2ir 



This compares fairly well with the value A\ ~ 0.25 GeV 



,||. Note that in Ref. |7[ the coefficient of A 1 /Q in Eq. ( pj|) 
6.3 rather than 7r 2 /2 ~ 4.9. This is because the calculation was performed there 
using the massive- gluon phase space with a massless matrix element, while here we include a gluon 
mass throughout. As discussed above, the fact that the correction to <tl is controlled by gluon 
fragmentation makes the calculation of its coefficient by the massive gluon technique less reliable 
than that of quark-dominated quantities. Nevertheless the technique is useful in revealing the 
existence of the 1/Q correction and its cancellation in cr to t, and it is reassuring that the calculations 
with and without mass effects in the matrix element give numerically similar results. 
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